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^ ■ Abstract 

o 

We investigate the propagation of a suddenly applied tension along a ther- 
mally excited semi-flexible polymer using analytical approximations, scaling 

j>- ■ 
00 

' arguments and numerical simulation. This problem is inherently non-linear. 

o ' 

We find sub-diffusive propagation with a dynamical exponent of 1/4. By gen 
■ eralizing the internal elasticity, we show that tense strings exhibit qualitatively 

different tension profiles and propagation with an exponent of 1/2. 
Characteristic for soft matter systems such as polymers or membranes is the often subtle 
interplay between energy and entropy [|]|§. Thermal motion determines the conformations 
of these systems crucially. For instance, the typical end-to-end distance even of a semi- 



o 



^ ■ flexible polymer (let alone a Gaussian polymer) is much smaller than its contour length. 
S-T 



Likewise for vesicles, thermal fluctuations can store a significant part of the "true" surface 
area. From the perspective of the projected length or area, this part is hidden. For vesicles or 
membranes, the hidden area can be pulled out without stretching the true area by application 
of a localized force such as the suction pressure in a micropipet [[| or the action of optical 
tweezers ||. For polymers, recent advances in experimental techniques using magnetic 
beads, optical tweezers or flow fields made feasible detailed experimental studies of the 
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conformations of single long molecules |5|-[T0|). With these techniques it has become possible 
to measure the force-extension curve for these soft systems and to compare them with 



theoretical predictions fTT|-|l4|| . These studies focus on the equilibrium conformations under 
the action of a stationary force. 

The purpose of this letter is to investigate the dynamics of the straightening of initial 
thermal fluctuations. The approach to the stationary state takes time because these objects 
are immersed in a viscous medium (usually water). Since we are interested in the principle 
mechanism, we will simplify the system as much as possible still retaining what we believe to 
be the salient features. Rather than studying the complicated dynamics of a two-dimensional 
membrane in three dimensional space, we will, therefore, study the straightening for a one- 
dimensional semiflexible chain confined to a two dimensional world. The unwinding of a 
Gaussian polymer in a homogeneous flow is the only other system so far for which the 



dynamical approach to the new equilibrium has been studied using scaling arguments |T2 
In contrast to this work, the conformation of our chains are dominated by bending modes. 
Also, we will consider the case of pulling one end holding the other end fixed, unlike the free- 
end situation in a homogeneous flow considered in [ |TT1 , |i"2"| . Our problem is thus analogous 
to the case of pulling on a vesicle containing a fixed volume of fluid [[|. 

Our main observations are that tension spreads slowly in a long chain, with a scaling 
law depending on the elastic character. We give numerical results along with simple scaling 
arguments to explain them. 

Model: To make the calculation tractable we will neglect the thermal fluctuations in 
our dynamics: their only role is to prepare a rough initial state at time t = 0, in which the 
total contour length per unit of projected real-space length is greater than one. We also 
assume that at t > the bending stiffness k is negligible compared to the applied force f; 
the only role of k is then to set the initial roughness of the contour (see below). 

Let us write the shape of a "polymer" confined to 2d as r(s) = (x(s),y(s)) along its arc 
length < s < L. We fix the right end to r(L) = 0. At the left end, a force f = — fit pulls 
in the — x direction. 



The total energy is then given by 

F = -fr(L) + l f L d S1 (s)r(s) 2 (1) 
I Jo 

where dots denote derivatives with respect to arc length. The Lagrange multiplier 'y(s), 
which can be interpreted as the local tension, is necessary to preserve the arc length con- 
straint i"(s) 2 = 1 which arises from the inextensibility of the polymer. We assume local 
isotropic dissipation and ignore long range hydrodynamic interactions. The equation of 
motion is then 

r-^r = -5F/5r = d s ( 7 r) + (f + jr)5(s) (2) 

where T is an inverse friction coefficient. 

There is no explicit equation of motion for the Lagrange multiplier 7. This quantity has 
to be determined from 

(l/2r)^r 2 = rd 2 ( 7 r) = 7 - 7 t 2 = 0. (3) 

At any time t, this relation for 7 is a second order differential equation which requires 
for its solution two boundary conditions. The first, at s = 0, is 

7(0) = / (4) 

The fixed boundary condition at s = L requires (l/T)<9 t r(L) = j(L)r(L) + r y(L)r(L) = 0. 
This relation implies both 

i{L) = (5) 

and r(L) = because rr = 0. The tension profile 7(5) can thus be calculated for any 
instantaneous conformation. This profile can then be put into the equation of motion for r 
which can then be integrated one time step. 

Linear theory: Even for our simplified model, the coupled non-linear equations (^j and 
(§) cannot be solved analytically. Let us consider first the linearized problem, which will fail 
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in an instructive way. Choose a single mode initial conformation of the form y(s,t = 0) = 
aP q sings. Note that x(s) follows once we specify y(s) because of the arclength constraint. 

Ignoring the non-linear term, the equation for the tension becomes 7 = which leads 
with the boundary conditions to the flat profile -y(s) = f. Thus the tension spreads in- 
stantaneously in this approximation. In the linearized equation of motion for y, the modes 
decouple and one obtains the usual relaxation form, a g (t) = a q > e~ T f q2t . 

The failure of the linear theory can be seen clearly by looking at the first correction 
term. Writing 7(s) = / + 7i(s), we have 71 (s) = fr 2 , which implies for the single mode 
configuration 71 (L) ~ — L 2 /|a°| 2 g 4 . This is less than the lowest order term / only if 

L<£l/q 2 \a° q \. (6) 

Thus depending on the wave-length and the amplitude of the initial conformation, there is a 
maximal chain length over which the linearized theory can be applied. For any longer chain, 
the problem becomes inherently non-linear. 

Simulations: To gain insight into the non-linear problem, we simulated a discretized 
version of the equations of motion. The initial conformation is written as 

2/0,0) = J2 a Q sin( l s ( 7 ) 

1 

with q = utc/L, n > 1. This choice guarantees that y(0, 0) = y(L, 0) = 0. For simplicity, we 
also restrict y(0,t) = for all times so that the left endpoint moves only horizontally [|i~5|j . 
The small initial amplitudes are chosen from a Gaussian distribution with width 

(\a° q f) = T/nLq\ (8) 

where T is temperature in units of energy and k is the bending stiffness. 

In Fig.l, we show ten snapshots of the configuration and the corresponding tension pro- 
files. The tension spreads inward, decreasing the amplitudes in this range. For a quantitative 
analysis, we define a penetration length £ (t) over which the tension has already spread. For 
concreteness we will define £ as the point where 7(£) = f/2. The initial value and the time 
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evolution of this penetration length depends on the initial conformation. Averaging 50 runs, 
we find the power law 

m) - t z (9) 

where the dynamic exponent z = 0.24 ± 0.01. We will now show how this empirical scaling 
law follows from a simple physical picture. 

Scaling argument I: An exact solution of the non-linear equations of motions looks 
impossible. Good insight, however, can be gained from the structure of the equation for the 
tension profile @. Since 7 has the same sign as 7, the two boundary conditions and @ 
imply that 7(3) is positive, monotonically decreasing and convex ]l(J. We will now seek an 
approximate effective equation for the time development of 7. For uniform tension the modes 
decay as a q {t) = exp(—Tjq 2 t)a q . Since this says that the short-wavelength modes decay the 
fastest, we will make the "adiabatic" Ansatz that 7 varies slowly in space compared to the 
wavelengths 2ir/q of the relevant modes; thus we write 

y( s ,t) = Y, a <i( s i t ) sin< i s (10) 

q 

and regard a q as slowly-varying functions of arclength s and of time: 

a q (s, t) = exp(-r7(s, t)qH)a° (11) 



Using eq. (]|), the average time-dependent local curvature then can be approximated as [17 



(y(s,t) 2 ) =£-^g 4 exp(-2r 7 M)g 2 t). (12) 

q KL, Q 

Written in this form, the curvature < y(s,t) 2 > depends only implicitly on the arc length s 
via the function j(s,t). With this simplification, Eq. (|3|) for the unknown time dependent 
tension acquires the simple form 

7(«) =7E:J^ 4 exp(-2r7(s)g 2 t) es -^M. (13) 

The equation for the tension profile has thus become a simple mechanical equation for 
the one-dimensional motion of a particle with unit mass in a potential ^(7) where 7 plays 



the role of position and the original spatial variable s plays the role of time. The original 
time t becomes a mere parameter. The potential exhibits two regimes, 



I 2 , 7 < 7c 

V(7) ~ { T (14) 

_ 7 3/2 ? 7>7c 



/t(n) 1 / 

separated by a cross-over tension 7 C = 1 /Tg^t. Here, g m is a high momentum cut-off for the 
g-modes. 

The particle starts with "initial" condition 7(0) = / and has to reach zero velocity at 
"time" s = L because of 'j(L) = 0. In the limit L — > 00, it then follows that the particle 
has zero total energy from which we obtain the initial condition 7(0) = — 2(V(/)) 1 / 2 < 0. 
Identifying 7(0) with — we obtain 

^(t) ~ {K/Tf/\Tftf^ (15) 

in good agreement with our numerical simulation. 

For a rough estimate of the relevant scales, let us ignore all factors of order unity. For the 
friction coefficient V we can take the inverse of the bulk viscosity 77, so V = 1/rj = 100cm 3 /erg 
sec. A typical weak force is f — 0.1 pN 0. For DNA, k/T ~ lOOnm |J and thus £(£) ~ 
l/im(t/sec) 1 / 4 . For actin, k/T ~ lO/xm |^8| , so £(£) ~ 10/im(t/sec) 1//4 . Thus, the dynamics 
of straightening in biopolymers should be accessible to video microscopy techniques. 

Generalization: In the analysis above we took the initial chain configuration to be 
governed by bending modes. Another case of interest is when a floppy chain is initially 
under tension, and the tension is suddenly increased at time zero. To cover both cases let 
us replace the sample for the initial amplitudes in Eq. (§) by (|ag| 2 ) = A 3 ~ 2b /Lq 2b , where A 
has the dimension of a length. Then the case 6 = 2, where A = k/T, was discussed above, 
while 6=1 (where A = T/E) corresponds to a string with tension E. 

Our scaling argument can be repeated for general b > 1/2. One finds the effective 
potential 

7(7) ~ _^3-26 (r/t)(2b -5)/2y,-l/2 (16) 



for 7 3> 7 C . This potential leads to a penetration length ~ t( 5 - 26 )/ 4 . 

We have studied the case 6=1 numerically. Here we find a dynamical exponent z = 
0.50 ± 0.01 clearly different from the prediction z = 3/4 which would follow from the above 
scaling argument. 

What went wrong? Closer inspection of the tension profiles indicates a qualitative dif- 
ference between the tense string b = 1 and the semi-flexible 6 = 2 cases. For the string, the 
10 snapshots shown in Fig. 2 and the corresponding tension profiles reveal that the tension 
profile is not exponential but rather decays quite linearly over a long region, and the moving 
boundary of this region remains sharply defined. Intuitively the difference stems from the 
fact that more of the initial excess length is in short- wavelength modes in the tense string 
case. Since these modes are damped the fastest, the string straightens immediately when it 
feels the tension, and so most of the resistance to straightening comes from pulling a straight 
string through a viscous medium, which gives a linear tension profile. The straight region 
ends at a point controlled by the total length pulled so far. In the semirigid case, most of 
the excess length is initially in long-wavelength modes, which decay slowly. The tension 
propagates forward before the string has a chance to straighten, and so its front is not so 
well-defined. 

Scaling argument II: We can turn these words into another simple scaling argument. 
We separate the string into two parts. In the left, straightened, part of length £(£), we 
assume that all fluctuations have already been pulled out. Then the tension profile becomes 
l( s ) — f~i(L)s. In the right, unperturbed, part, pulling has not yet had a significant effect. 
Thus the tension profile for s > £(t) is exponential, 7(5) = 7(£(t)) exp(— (s— with 
a time-independent penetration length £ init = (1/r 2 ) 1 / 2 as given by the initial configuration. 
Matching the two profiles so that 7 is continuous at s — £(t), we obtain 

7 (£(f)) = //(I + m/iinit) ^ Kinit/m- (17) 

The rate d£(t)/dt with which the straightened regime grows is proportional to the force 
7 (£(£)) with which this regime "pulls" at the left end of the unperturbed regime. We thus 
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obtain d£(t)/dt ~ 7 (£(£)) ~ f£init/£(t). This equation is readily solved to give 

C(t) ~ (ft) 1 / 2 . (18) 

This "diffusive" behavior indeed corresponds to our simulation results for 6 = 1. We can 
easily understand it in terms of the intuitive picture sketched above: the total friction on 
the straight segment is proportional to its length £(£), and so the string velocity is v ~ 
Letting a be the initial excess contour length divided by L, conservation of string says that 
the front velocity is £ = v/a ~ Solving this equation then reproduces £ ~ \ft. 

In Fig. 3, we show the numerically determined exponents for various values of 6. We have 
strong numerical evidence that z — 1/2 for 6 < 1 and z < 1/2 for 6 > 1. This plot shows 
that the two scaling theories yield upper limits to this exponent. 

Conclusion: We have seen how thermal fluctuations, together with a viscous surrounding 
medium, impede the transmission of a suddenly imposed tension in a polymer. As we 
pointed out, this phenomenon is inherently nonlinear. In a model neglecting stochastic 
noise, simulations show that the tension propagation in a stiff chain is subdiffusive with an 
exponent 1/4. We showed how this law also follows from a simple self-consistent scaling 
theory. We also generalized both the model and the scaling theory to the case of chains with 
arbitrary initial elasticity. For tense strings, the spreading becomes diffusive, in the sense 
that the tension is nonnegligible in a region of size ~ y/i. The application of these ideas to 
the case of suddenly-pulled membranes || will be addressed in future work. 
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This work was supported in part by the US/Israeli Binational Foundation grant 94-00190 
and NSF grant DMR95-07366. 
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FIGURES 

FIG. 1. (top) 10 subsequent snapshots of conformations (x(s),y(s)) of a semi-flexible chain 
(b = 2) with total chain length 250. The snapshots are taken at equidistant times. The inset shows 
how two adjacent bends are straightened, (bottom) Corresponding tension profiles as a function 
of the arc length. 

FIG. 2. Same as Fig.l for a string (6=1). 

FIG. 3. Dynamical exponent z as a function of b. Error bars have the size of the symbols. The 
lines refer to the two scaling estimates. 
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Fig. 1 (Seifert et al.) 
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Fig. 2 (Seifert et al.) 
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